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Abstract 

We prove that for quantum lattice systems in d < 2 dimensions the addition of quenched disorder 
rounds any first order phase transition in the corresponding conjugate order parameter, both at 
positive temperatures and at T = 0. For systems with continuous symmetry the statement extends 
up to d < 4 dimensions. This establishes for quantum systems the existence of the Imry-Ma 
phenomenon which for classical systems was proven by Aizenman and Wehr. The extension of the 
proof to quantum systems is achieved by carrying out the analysis at the level of thermodynamic 
quantities rather than equilibrium states. 
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I. INTRODUCTION 



Quenched disorder is known to have a pronounced effect on phase transitions in low 
dimensional systems. As Imry and Ma pointed out in 1975 jl|, disorder can prevent the 
appearance of discontinuities, and of long range order, associated with first order transitions. 
Nonetheless, for a number of years there was considerable uncertainty about the mechanisms 
involved in such a "rounding effect" , and consequently also of the circumstances in which 
this effect appears j^J. Rigorous results in a series of works in the 1980s {3, 4] established the 
conditions for rounding in classical systems. The situation in quantum systems remained 
uncertain, with recent suggestions jf]] that for such the effect might not have the same general 
character as in the context of classical statistical mechanics. 

In this work we present a modified version of the argument of Ref. |4j which allows us 
(as previously announced in Ref. [^) to extend the classical results to quantum systems. It 
was not obvious that this should be possible, since the statistical mechanical description 
of a d dimensional quantum systems often has the appearance of a {d + 1) dimensional 
classical system, and in addition some of the tools used in the classical case have no clear 
quantum generalizations. The latter difficulty is resolved here by focusing on thermodynamic 
quantities, bypassing some of the subtle issues related to equilibrium states which played a 
role in Ref. |4|. 

Before presenting the formalism in which the general result is expressed, we start with 
some specific examples to which the results apply. 



II. EXAMPLES OF THE ROUNDING EFFECT 



A. Transverse Field Ising Model 

The transverse-field Ising model in a random longitudinal field is defined by the Hamil- 
tonian 

^tfim = _ j a3 xa3 y <7i,x - J> + ^W*, (2-1) 

&,y) 

where Pauli matrices describing the j = 1, 2, 3 components of a spin at site x, and 

the first summation is over pairs of neighboring sites x, y G Z d . The symbol r\ x represents 
a random field in the j = 3 direction, whose values at different sites are given by indepen- 
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dent and identically distributed (i.i.d.) random variables, satisfying certain mild regularity 
conditions (to be specified in Section HV Bl below) . 

In d > 2 dimensions the nonrandom version of the system (with e = 0) exhibits spon- 
taneous (or residual) magnetization at h = if J/A and J (3 are large enough. In d = 1 
dimension this occurs only at zero temperature, T = = 0, i.e. in the ground state. More 
explicitly, existence of spontaneous magnetization means that 

M{0) := ^[Mo+(/3) - M -(/?)] > (2.2) 

where Mq±{0) are the two limits of the Gibbs state expectation values in finite (rectangular) 
domains T C Z d , 

M 0± (/3):= Km lim -- ( . (2.3) 

Since the h limit is taken last, the boundary conditions on the finite systems do not affect 
the value of the limits in ( I2.3P (for the convenience of translation covariance our default 
choice is periodic boundary conditions). In view of the symmetry of the model, M 0+ = 
-M _ = M(/3). 

The general result which is proven below (Theorem HV.lj) implies: 
In d < 2, the system described above, with the Hamiltonian ( 12. ip . at e ^ has M((3) = 
at all (3 < oo. In particular, even if there is more than one ground state, they all yield the 
same bulk-average value for the mean magnetization. 

This can also be restated in thermodynamic terms: Letting ^(h, 0) denote the ther- 
modynamic limit of the specific free energy, which is defined below in ( 13. 6p . the different 
magnetizations can be expressed as directional derivatives [3] 



Mo±C9) 



dh± 



(2.4) 

h=0 



Non-vanishing spontaneous magnetization corresponds therefore to a discontinuity of a first 
order derivative of the free energy, here at h = 0, and thus to a first order phase transition 
in the terminology of Ehrenfest [8]. 

Note that the order parameter for the phase transition is the expectation of the volume 
average of the local operator a^ yX which in the Hamiltonian (12. ip appears coupled to the 
random field. Such conjugacy is the essential condition for the general result presented here. 
The theorem does not imply that the rounding of the magnetization in the '^-direction' 



would occur when random terms are added only to the transverse field A. Indeed, as has 

n 

long been known[9|, when only such disorder is present spontaneous magnetization does 
persist in two dimensions. Such transversal disorder can have other subtle effects on the 



ferromagnetic-paramagnetic transition 
work. 



10], but these are beyond the scope of the present 



B. Isotropic Heisenberg model 

The isotropic, nearest-neighbor Heisenberg ferromagnet in a random magnetic field is 
described by the Hamiltonian 

•H Hcis = -jJ2 & *- ff y-J2& + tf*)-Vx (2-5) 

where J is a positive real number, a x is a vector spin operator associated with the site x G Z d , 
and ff x is a magnetic field which varies randomly from site to site. We now assume that these 
random field variables are not only independently and identically distributed, but also that 
their distribution is rotation invariant. We also assume that ff x is nonzero with probability 1. 



In the case of a uniform magnetic field (e = 0) it is universally believed that for d > 3 
this system exhibits spontaneous magnetization at h — below some J-dependent critical 
temperature (although a rigorous proof is lacking). Letting (-) r now denote the Gibbs 
state of this system on a finite domain T C Z d with periodic boundary conditions at inverse 
temperature /3, spontaneous magnetization can be expressed, similarly to the previous cases, 
by the statement that 

M„(/3) := Urn to ± ^0 (2.6) 



where e is any fixed unit vector (indicating the direction along which h is taken to 0). 
In terms of t 



In terms of the thermodynamic limit of the specific free energy the same quantity can be 

m 



Note that as in the previous example the randomness is conjugate to the quantity whose 
density is the magnetization considered here. 



Theorem IIV.2t below, has the following implication for this system at e ^ 0: 
In d < A, the random-field Heisenberg ferromagnet described above has Mq(/3) 
hence J-(h, /3) is differentiable in h at h = 0) for any (3 including (3 = oo. 



(and 



As will be seen in Section IIVCI the vanishing of the spontaneous magnetization implies 
also the absence of ferromagnetic long range order. 

III. A GENERAL FORMULATION 
A. The system and its Hamiltonian 

We consider here systems on homogeneous d dimensional lattices, which for simplicity 
we take to be Z d . Associated with each lattice site x G Z d is a quantum system whose state 
space is isomorphic to a common finite dimensional Hilbert space Hq. The local systems 
are coupled through a Hamiltonian which for a finite region T C Z d , with free boundary 
conditions, takes the form 



Here Qa is an operator which acts on the quantum degrees of freedom in A (i.e. it is 
described by an operator acting in the space ® x &aH x ). The interaction is assumed to 
be translation invariant (Qt x a = U^QaU x , with T x denoting translations on the lattice 
and U x the corresponding unitary operators). Unless stated otherwise it is assumed here 
that the interactions are of finite range (i.e., Qa = for all sets with diamA > R, at 
some finite R). However the stated results apply also to a a broader class of long range 
interactions with a sufficiently fast power law decay, which is described in Appendix |A] In 
particular, Theorem II V. 1 1 applies to two-body interactions satisfying ||Q x ,j/|| < C/\x — y\ 3d ^ 2 , 
and Theorem IIV.2I applies if ||Qx,y|| — C/\x — y\ a with some a > d — 2. 

Disorder is expressed in (13.11) through the random coefficients r) x . Their variance is fixed 
at 1, so that the strength of the disorder is controlled by the parameter e. The terms 7] x 
multiply are assumed to be of the form k x = U^.kqU x , with kq an operator which acts on a 
finite cluster of sites, whose size may be greater than 1. It should be noted that {t] x } appear 
in (13.11) as additions to a uniform parameter h. This parameter plays an essential role for 
the results presented here. The symbol r] without the subscript indicates the collection of 




(3.1) 
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random fields at all sites, and rfp denotes the random fields in some finite region T C 7L d . In 
systems with several families of independent random terms our results hold for the observable 
associated with each family considered separately, but we will not discuss this point here 

(see Q). 

Hamiltonians with other boundary conditions, H^g V , are defined analogously, so that the 
terms within T are the same as for free boundary conditions, while at a finite distance from 
the boundary there can be additional terms which are uniformly bounded in norm and in 
the size of the clusters of the directly affected sites. Among the allowed options are periodic 
boundary conditions, denoted by£> = per, which will be the default choice if no other is 
explicitly specified. 

As usual, the finite volume partition function is denoted 

Z r>B (h,f3,e;r)) := Tr exp(-/?i^f), (3.2) 
and the corresponding Gibbs state is 

Tr f ... e-^r'A 
Zr, B {h,/3,e;r)) 

When the values of the superscripts is deemed obvious from the context, they will sometimes 
be omitted. 

B. The quenched free energy and its infinite volume limit 

The free energy for a finite system is denoted here by 

F TtB (h, /?, e; rf) := y log Z r (h, /?, e; 77) , (3.4) 

We will mostly be concerned with square domains = [-K, K] d D Z d with the periodic 
boundary conditions, and so we introduce the abbreviation 

Fr{v) = F rK>pei (h,(3,e;r)). (3.5) 
The free energy density, per unit volume, is 

Fr, B (h, P, e; V) ■= Ft Ah, P, e, v)/\T\ • (3.6) 

where |T| is the number of sites in T. As the notation suggests, for finite systems the free 
energy depends on the choice of boundary conditions and on the disorder variables. In the 
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thermodynamic limit, however, the dependence of its density on the boundary conditions 
disappears and so does the dependence of its typical value on the disorder variables. The 
following generally known result (we refer here to such statements as Propositions) applies 
to all systems of the type described above. 

Proposition III.l. If the random fields are of finite variance and form a translation in- 
variant and ergodic process, then for any ft G [0, oo] there is a full measure set M of field 
configurations for which the infinite volume limit 

T(h,ft,e):= Km JF Tl ,b{K ft,e,rj) (3.7) 

L— >oo 

exists for all h and its value is independent ofrj and the boundary conditions B. 



Since much of the early discussion of this result [12] was limited to the classical case let 
us add that the argument used in that context extends directly also to quantum statistical 
mechanics. In essence: the free energy density can be approximated up to an arbitrarily 
small correction, of order 0(Lq x ) (see Inequality (16. ip ) by that of a system obtained by 
partitioning the space into large blocks (of length L Q ) among which the couplings were 
removed. For the approximants the density converges by either the law of large numbers, in 
case of independent {t] x }, or by the ergodic theorem in the more general case stated above. 
That implies Proposition IIH. l] 

The following remarks are of relevance: 

i. Since convexity is automatically inherited by the limit, Fih^ft.e) is concave as a 
function of h. 

ii. The case ft = oo corresponds to the ground state energy: for almost every rj the limit 

J*(/i,oo,e) := lim F(h, ft, e) (3.8) 

exists and is equal to the limit of the energy densities of the finite volume ground states of 
the random Hamiltonian (whose value is independent of the boundary conditions). In other 
words, for the free energy the limits ft — > oo and L — > oo are interchangeable. 

Hi. The uniqueness of the free energy density does not extend to uniqueness of the Gibbs 
states, or ground states in case ft — oo. The question which our results address is whether 
the different Gibbs states (or ground states) of the given Hamiltonian can differ in their 
mean magnetization, that is in the volume averages of (k x ). 



iv. The following condition will allow us to weaken other assumptions. We will say that 
our system satisfies the weak FKG condition with respect to k if 

(^>? (V) > (V') (3.9) 

whenever x G T, h > h', and r\ y > n' y for all y G T. This condition is known to hold in the 
transverse field Ising model described in Section Til Al {9]. 



C. Quenched disorder with continuous symmetry 

For systems with a continuous symmetry, the Hamiltonian (with the free boundary con- 
ditions) is of the form 

Act xgf 

where the Q interaction is not only translation invariant but also rotation invariant in the 
following sense. We assume that the group of rotations SO(N) has a representation by local 
unitary operators whereby for each R G SO(N) there is a corresponding R x acting in r H. x . 
The interaction Q then is rotation invariant if and only if 

qa= (n^Wfn^) ( 3 - n ) 

for all finite A C Z d and all R G SO(N). At each site x G Z d instead of a single operator 
k x there is now a collection of operators k x = (k x ,i, k x ,n) which under the above action 
of SO(N) transform as a vector. That is: for each R G SO(N) and v G M> N 

R-\v-k x )R x = (Rv) ■ H x . (3.12) 

Hamiltonians with other boundary conditions, Gibbs states, etc. are defined in analogy 
with the previous section. (In the classical case 4J it was also necessary to assume that 
the interactions transformed smoothly under nonuniform rotations, however for a finite- 
dimensional local state space this is automatically true.) 

Likewise, the random fields are now given by random vectors ff x in M. N . Clearly rotation 
invariance is lost in the presence of such terms in the Hamiltonian. However, we will assume 
that the symmetry is retained in the distributional sense, i.e., that for any rotation R the 
random variables ff x and Rff x have the same distribution. 
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IV. STATEMENT OF THE MAIN RESULTS 



A. Two perspectives on 1 st order phase transitions 

As was done in the context of the examples of Section [Til the general statements will be 
presented in two equivalent ways: in their thermodynamic formulations, which is the level 
at which the results are derived here, and then in the language of statistical mechanics, 
i.e. properties of Gibbs states. Since the results concern (the absence of) first order phase 
transitions let us recall the latter's dual manifestations. 

1. In thermodynamic terms, a 1 st order phase transition is associated with the disconti- 
nuity of the first derivative of the free energy with respect to one its parameters. By 
default this parameter will be denoted here by h. 

2. In the terms of statistical mechanics a 1 st order phase transition is expressed in the non- 
uniqueness, among the infinite volume Gibbs equilibrium states, of the bulk density of 
some extensive quantity. Our discussion concerns the case when this is the quantity 
whose coupling parameter in H is the field h which is randomized by the disorder. 

One occasionally finds that the equivalence of these two statements is not fully trusted in 



the present context (of disordered systems) 13[]. We shall therefore briefly recall below, in 
Section TlV CI some pertinent known results. 



B. Thermodynamic formulation 

Our first result concerning systems of the type described in Section IIII A\ with rj x given 
by independent and identically distributed (i.i.d.) random variables, is: 

Theorem IV. 1. In dimensions d < 2, assuming the variables 7] x are i.i.d. with absolutely 
continuous distribution and with more than two finite moments, the quenched free energy 
density J 7 is differentiable in h at all values of h, e ^ 0, and (5 < oo. 

Furthermore, for /3 < oo, the assumption of absolute continuity can be relaxed, requiring 
instead that the distribution off] has no isolated point masses, or alternatively that the system 
satisfies the weak FKG property with respect to k. 
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Let us recall that a probability distribution of a variable q is said to be absolutely continu- 
ous if it is of the form p(q)dq, with a density function p(q). In particular, the corresponding 
measure on R has neither 'delta function' terms (point masses), nor any component which 
is supported on a Cantor type fractal set. Point masses are isolated if they are separated 
from the continuous part of the distribution. The statement that the random variables rj x 
have more than two finite moments means that 

Av|7fe| 2+a < oo (4.1) 

for some 5 > 0, where Av denotes the average over the random fields. 

As in the classical case, in the presence of continuous symmetry, as understood above, 
the 'rounding effect' extends to higher dimensions: 

Theorem IV. 2. In dimensions d < A, any isotropic system of the type described in Sec- 
tion \IIIC\ with the random terms being independent with an identical, rotation invariant 
distribution with more than two finite moments and Prob(r/ = 0) =0, the quenched free 
energy density T is differentiable in h at h = 0, for any e ^ and (5 < oo. 

For (3 < oo, the conclusion still holds if there is a nonzero probability of a random term 
being zero. 

We note that an apparently weaker condition on the distribution of ff compared to The- 
orem IIV.1I is adequate because what will ultimately be important is the distribution of the 
component in an arbitrary direction. With the assumption of an isotropic distribution for 
the vector, the components can easily be seen to satisfy the stronger conditions used in 
Theorem IIV.1I or Proposition IVII.2I 

C. Statistical mechanical implications (no long range order) 

In this section, we establish the relationship between phase transitions (understood in 
terms of non-differentiability of free energy) and long range order, culminating in corollaries 
of Theorems IIV.1I and IIV.2I which reframe the results in statistical mechanics terms. 

We say that a system exhibits long range order with order parameter k (for some set of pa- 
rameters and disorder variables) if the mean value of the bulk averages k\ = |A|~ X XLga k xi 
for A — > Z d , depends on the boundary conditions B with which the infinite volume state is 

10 



constructed. That is: if there are two sequences of cubic domains A and T increasing to Z d , 
and two sets of boundary conditions which yield different values for 

lim lim (R,\) r R . (4.2) 

When the sensitivity to the boundary conditions is found for the full bulk averages, i.e. 
if different values occur for 

Km (K V ) r>B , (4.3) 
we say that there is long - long range order. 

Remark: If for a given system there are multiple limits of the expectation values in (14.21) . 
then the system has multiple infinite volume KMS states [l^j]. Convex combinations of the 
different states (which will also form KMS states) will exhibit non-decaying correlations of 
k: {k x k v ) — (K, x )(n y ) -/> as \x — y\ — > oo. The latter condition provides another aspect of 
("short"-) long range order, which is often used as its definition. 

The basis of the connection to the thermodynamic quantities discussed in the previous 
section is the relation 

— dh — = ifiZ^^/r.flW- ( 4 - 4 ) 

A useful implication of convexity is that if a sequence of differentiable convex functions, 
such as T of Theorem IIII.ll converges pointwise (in h) then also their derivatives converge 
to the derivative of the limiting function, wherever that function is differentiable. Without 
assuming differentiability of the limiting function one may still conclude that all the deriva- 
tives' accumulation points lie in the interval spanned by the left and right derivative of the 
limiting function. (An elementary proof of that can be obtained by considering the relations 
among quotients of the form [J 7 ^) — J r (hi)]/[Ii2 — ^-l] for suitably chosen collections of 
intervals [/ii,/^]-) This has the following relevant implication. 

Proposition IV. 3. Under the assumptions of Proposition \TTLl\ for any set of the parame- 
ters (/?, h, e) at which J 7 is differentiable in h 

li m _!_ V" U x )f R (77) = ^- (4.5) 
1 1 xer L 

for almost every realization of the disorder!], and any choice of the boundary conditions B. 
Furthermore, also: 

1 dJ- 

x€T L 
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Without assuming differentiability, one may still conclude that 

{ K x)r L ,B [V) e 

x£T L 



LIM -— V ( Kx )$ b 0?) e 



dh- ' d/i- 



(4.7) 



LIM denoting the collection of accumulation points, for different boundary conditions, 

d 
dh± 



and possibly different sequences of volumes, and -Jf-r denoting one-sided derivatives with 



respect to h. 

Of particular relevance for us is the conclusion that when J 7 is differentiable there is 
only one possible value for the limit (14.2jl . which by Proposition IIH. ll is independent of the 
boundary conditions. For classical systems more can be said: if J 7 is differentiable not only 
do the mean values of the observable Rr converge, but the distribution of this quantity with 
respect to the Gibbs state collapses onto a point. Such a stronger statement concerning 
quantum fluctuations is not known to be true. Nevertheless, Equation ( 14. 5 p holds also for 
the quantum systems. 

Proof. The relations ( 14. 5 p and ( 14. 7p follow directly by the above stated property of convex 
functions (i.e., the implication for the derivatives of convex functions' pointwise conver- 
gence) . 

To prove ( 14. 6ft and the corresponding extension of ( 14. 7ft . let Fp' B ' A denote the free energy 
with the fixed field within A changed by 5, so that 



FT ( k *)v k ,b (v) 



1 ur T K ,B 



(4.8) 



i r I ^ i,1 ^ w; Til 88 
A standard estimate, which is presented in Proposition IVI.ll below, implies that 

*r£? - F Lb = Fg£ - FLo + (W V r\ II) (4-9) 

uniformly in K, where V'P denotes the terms in the Hamiltonian crossing the boundary of 
T L , i.e. 

A: AnT <£{(/), A} x&doV 

Then for any rj in the full measure set M on which the limits in Theorem IIH. II are defined, 
this implies that 



lim - F* KtB ) I \V L \ = F{h + 6,0)- JF{K 0). (4-11) 

^oo ; K>L V ' / 
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each of these differences is a concave function of 5, so the above mentioned convexity argu- 
ment gives 

Pi . r FIT FlT 1 

(4.12) 



LIM 4^r M 'B /IIlI e 

L,K^oo;K>L 35 Vk ' B 



3h- 3h+_ 

and the desired results follow immediately from Equation (14. 81) . □ 

In particular, the results which are formulated in Theorems IIV.1I and IIV.2I as statements 
of differentiability of the free energy directly imply the following. 

Theorem IV. 4. In dimensions d < 2, under the assumptions of Theorem \IV. 1\ the system 
almost certainly does not exhibit (short or long) long range order. 

Furthermore, for systems with continuous symmetry which meet the assumptions of The- 
orem UVJ^ in dimensions d < 4 almost certainly there is no (short or long) long range order 
ath = 0. 



V. THE FREE-ENERGY-DIFFERENCE FUNCTIONAL 

The above discussion reduces our main results to Theorems IIV. II and IIV.21 A key tool for 
the proof of these theorems is a quantity which describes the differences in the finite volume 
content of free energies between two of the system's equilibrium states, which are constructed 
to have the two extremal mean magnetizations, and J^, assuming they are not the 
same. The difference, which is denoted below as Gl, will be shown to satisfy contradictory 
bounds: i. an absolute upper bound on \Gjj\ which is due to the observation that the two 
free energies cannot differ by more than the magnitude of the interaction across the finite 
volume boundary, ii. a lower bound on the fluctuations of this quantity which reflects the 
idea that with a systematic difference in the magnetization the states' free energies will 
respond differently to the fluctuating random fields. 

More specifically, the free energy differences (properly defined) will be shown to satisfy: 

\G L \ < CL^ + DL^ 2 , (5.1) 

with D = for finite range interactions. In the presence of continuous symmetry for finite 
range interactions the bound is improved to: 

\G L \ < CL d ~ 2 . (5.2) 
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For the opposite bound it will be shown that if there is a first order transition then Gl/ ' L d l 2 
converges to a normally distributed random variable with variance b > 0: 

' -+N(0,b). (5.3) 



L d/2 



The construction of a quantity with the above characteristics is the main subject of this 
section. The contradiction between the two bounds yields the main result. 

Estimates similar to the above form the basis of the argument of Imry and Ma[l|, and 
the more precise argument outlined above is similar to the one which was employed for 
the analysis of in the classical case ^|. However the technique introduced there for the 
construction of the auxiliary quantity Gl employed probability measures over states, defining 



what has since been called metastates 



15[. Unfortunately, although metastates can be 
defined for quantum systems [16| , not all the steps taken in Ref . |4j for the construction of 
Gl have such an extension. The progress made in this work is enabled by the observation 
that one can define a suitable quantity Gl in terms of just free energy differences, avoiding 
more delicate issues of quantum states. This makes it possible to formulate a proof parallel 
to the classical case (and indeed the nonrigorous arguments which inspired it). 

As a simple illustration of the concept one may first consider the zero temperature case of 
a system such as the random field Ising model with two distinct states, labeled by + and — , 
of different mean magnetizations and with only the 'canonical dependence' on 77 (cf. 4|). 
In that case, for Gl one may take the difference in the two states' energy contents in finite 
volumes: 



E L (v) ■= E+M - E- L (rj) = e ( ^ t} x k 3 
« 2e ( Vx j M ^2 eML d/2 N(0, 1) 




(5.4) 



va;er L 



where the last step, in which the central limit theorem is invoked, is valid only to the extent 
that one may ignore the adaptation of the system to the disorder. (A more precise statement 



is possible taking into account the dependence of the ± states on r\ {4]). The above provides 
only a suggestive example in the context of a special (and classical) system at T = 0. We 
turn now to a more general definition which would be suitable for our purpose. 

Our choice of Gl{v) f° r the general case is based on the following free energy 'second 
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difference' (in the field h, and disorder t]r L 



G 5 L , K M ■= 5 Av [F h K +5 { V ) - F^(r L ( V )) - Ft S (v) + F h K ~\r L { V ))\ VL ] 



(5.5) 



where rjj\ is the field obtained by setting r] x = for all x G T^, tjl is the restriction of 77 
to T/, and Av[-|r/i] is the conditional average over the random terms outside T^, i.e. with 
t]y l held constant. It should be noted that Fr (which is defined in (I3.5p ) refers to the free 
energy with the periodic boundary conditions. This choice assures translation covariance, 
which is mentioned below (condition 4. in Lemma ( 1V.1I) ) and used in the argument. 

When the double limit lim^o h m ir->-co G S L k (tjl), exists (for all tjl), the quantity defined 
by it has properties we desire of Gl{t]i). (The order of the limits is important here: taking 
8 \j after the infinite volume limit serves to construct the general version of the + and — 
states.) Inconveniently, the limits are not generally known (or expected in full generality) 
to exist. However compactness arguments can be applied to prove sufficient convergence 
along subsequences. The essential statement is the following. We use here the £ 1 -Lipschitz 
seminorm 1 1 ■ 1 1 (in lieu of a uniform bound on the derivative) which is defined as 

\m-nv')\ 



sup 

v,ri'es WV ~V Hi 

0< 1 1 77 — 77' 1 1 j < OO 



(5.6) 



with 



Lemma V.l. For any (5 < 00, there are sequences K~ — > 00 and 5i — > such that the limits 



lim lim f G 



T L,Kj 



Av 



% — yoo j— ¥00 

exist for all L and r\, and have the following properties: 



G 



<5 ! 

L,Ki 



(5.7) 



1. AvG L (r] L ) = 

2. Gl(til) depends on the values of r\ x only for x G 



3. IIIGrlll <e 



4. Av [G L \r] A ] = G V (T X 7]) whenever T„ X T L , = AcT L 



5. Gi (that is, Gl with L = 1, which is a function of one variable) has a distributional 



derivative G[ satisfying 



dGi 
diio 



< /3e 2 and 
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6. AvG' 1 M = f(^-S 

Furthermore, if the system satisfies the weak FKG condition with respect to k, then G\ is 
monotone. For f3 = oo the same results hold assuming that the distribution of i] is absolutely 
continuous. 

The proof is postponed to section IVIIIt as it is based on compactness arguments par- 
alleling those in Ref. |4|. Before turning to it, in the next two sections we will show that 
any such Gl{j]l) must obey the upper bound (15. ip . and where appropriate (15. 2p . and also 
satisfies (15. 3ft . 



VI. UPPER BOUNDS ON THE FREE ENERGY DIFFERENCE 

We will frequently use the following estimate to control the effects of changes in the 
Hamiltonian on the free energy. 

Proposition VI. 1 For any Hermitian matrices C,D of the same finite size, 

|logTre c - logTre D | < ||C-D|| (6.1) 

where || • || is the operator norm. 

As a simple application, this means that the change in free energy for the system defined 
by the Hamiltonian (13. ip when boundary conditions are added is bounded in terms of the 
total of the norms of all terms in the interaction which cross the boundary. 



A. A general surface bound 

The following lemma shows generally that the free energy difference due to random field 
fluctuations in a finite region can be bounded by the norm of the interaction terms crossing 
the surface of that region, and that for short range interactions this is proportional to the 
area of that surface. 



Lemma VI. 2. Assuming the interactions are either of finite range or satisfy flA.ip (of 



Appendix \A\) . for any pair of sequences Si — >■ ; Kj — > oo for which the following sequence 
of functions converges for all rjL 



G L ( VL ) = hm lim (G% K - Av \d% K ] ) (6.2) 
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the limiting function satisfies (with some C,D < oo): 

\G L \< CL d - x + DL d/2 . 

Proof. We will show that 

\G{ K \< cL^ 1 + bL d/2 + 0(5L° 



(6.3) 



(6.4) 



with c, b independent of L, K, and 5. The desired result holds then with C = 2c, D = 2b. 

To see this, we let be the smallest cube in Z d so that k x acts within for all x G Tl, 
and then split the Hamiltonian into contributions from Al and Tk\Al and terms connecting 
the two. We consider the quantity 



F£ |i (77):=--logTre X p(-/3 



:\lvd ■- —q lu s ^ ^ y~H [ h aZo + H Tk\a l ^ ) > ( 6 - 5 ) 
where the subscript refers to free boundary conditions, and the subscript * refers to periodic 
boundary conditions on the edge of Fx and free boundary conditions on the edge of A^,. The 
Hamiltonian in Equation (16.51) differs from only in the absence of terms connecting A^ 
to the rest of IV, which as in Equation (I4.10p we denote by . Then by Proposition I VI. ll 

(6.6) 



|^i^)-^)|<IKI 



and therefore 



< 



Ft 5 (v)+Ft 5 (r L (r,))\ 



7 h+5 
K\L 



fo) - F h K t{r L {r))) 



7 h-5 
K\L 



(V) + F^(r L ( V )) 



(6.7) 



Each Fk\l is a sum of a part from A^ and a part from the rest of Tl] the latter cancel 
each other exactly, and the former can be arranged in two pairs which differ only in the 
constant field term which has norm of order SL d . Taking the average conditioned on t]y l 
and noting that the definition of A^ makes Vjl independent of 7/r L , 



< 2Av||V^J 



0(5 L° 



'Ail 



for all rj, and since the 5 term is uniform in K, 

\G L ( VL )\<2Av\\V j 
This reduces Inequality (16. 3p to a comparable bound on , i.e 

AvIIVj II < cL^ 1 + bL d/2 . 



(6i 



(6.9) 



For finite range interactions it is obvious that Av 1 1 Vl 



(6.10) 

< cL d ~ l , and Appendix [X] 



discusses more general conditions under which Inequality (16.101) holds. 



□ 
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B. An improved bound for systems with continuous symmetry 



For systems with continuous symmetry the stronger bound of Inequality (15. 2p can be 
obtained through the analysis of soft deformations. The argument is similar to the one 
which was carried out in the classical context in Ref. |4j. Although the estimate may at first 
sight appear to be simply a calculation of the Bloch spin wave energy, that alone would not 
yield the desired result since the first order term is potentially much larger than the required 
bound. However, by comparing the distortion energies of opposite deformations one finds 
that in any situation for one of the distortions the first order term is of the desired sign. 
Therefore the analysis can be limited to the second order term - and that yields the claimed 
bound. 

We need to modify the definition of Gl slightly to accommodate random vector fields. 
The argument will be based on the effect of the fluctuations of a single component, described 
by expressions where the others have been averaged out. We let 

G%k(Vl) = \ Av [Fg(fft - F s K \r L {ff)) ~ F^(v) + F^(r L (ff))\e ■ fj L ] (6.11) 
G L ( VL ) = lim lim (G 5 £% (rj) - Av \d 6 £% .(e ■ ff}} ) (6.12) 

i— too j — too \ J L J J / 

where e is an arbitrary unit vector; note that the rotation symmetry of the system and of 
the distribution of ff mean that the right hand sight of the last expression is independent of 
e. 

Lemma VI. 3. For any system with continuous symmetry and isotropic disorder as described 
in Sectioning of interactions which are either of finite range or satisfy ( 1A.7j) . at h = 0: 



\G L {r]L)\<CL d ~\ (6.13) 
with C < oo, uniformly in L and tjl- 
Proof. We shall derive (16. 131) through uniform bounds on 

9 S L %(V) := Av [F$(rf) - F^(ff)\e ■ ff L ] , (6.14) 
for which we shall show that Assumption (1A.7I) implies 

\9% K m =0(L d ~ 2 ). (6.15) 
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The two quantities are related through 

1 



G L (rj L ) = ^ lim hm (g^ K .(v) - Av^i (e • 77) j , (6.16) 

and hence (16.151) implies the desired statement. 

Let p be the generator (in so(N)) of a rotation in a plane containing e, and for each 
x G Z d let p x be the generator of the corresponding rotation in the single-site algebra A x 
(Unlike in Ref. If], we will use the "mathematician's" convention that rotations are given by 
e ep , so that p is an antihermitian matrix and p x is an antihermitian operator). We introduce 
the slowly varying angles 

0, x e v L 
Qr := \ Wpi^ o < < L ) 

7r, ^l(x) > L 



(6.17) 



where di{x) is the distance from x to in the largest- component metric, i.e. 

d L (x) := min \\x - yW^ . 
We also introduce the associated rotations on fields and on Ar defined by 

R x ■= e 9 * p 
(Re(v)) x = RxVx 



Ro 



"xPx 



(6.18) 

(6.19) 
(6.20) 
(6.21) 



where Y is the relevant finite subset of Z d which can be inferred from the context. 
Re is unitary, and so we can rewrite the free energy F^ 5e (ff) appearing in f!6.14p as 



FkHv) = -^logTrexp {-^H^R, 



At this stage we could carry out a spin-wave analysis, writing 



Av [F^{ff)\e-fl L ] = --Av [logTrexp (-/3[Hff + AH e 



which by Lemma [VI. II implies 



gfM)\ < \\AH e \\ . 



(6.22) 



(6.23) 



(6.24) 



However this will not lead to the desired result, since HAif^H includes 'first order' terms 
(in V#) which scale as L d ~ x as well as the 'second order' terms (and higher) which scale as 
L d ~ 2 . To cancel the former we rewrite F^ Se (ff) using an opposite rotation and average the 
two expressions, to obtain 
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Av [F^(ff)\e-ff L ] = -^Av [logTrexp (-/3[Hff + AH B ] 

+ log Tr exp (-(3[H^ + AH 



(6.25) 



A combination of the Cauchy-Schwarz inequality, the Golden-Thompson inequality 
and Lemma I VI. II yields the general relation 



l0g Tr ^ ~ \ bg Tr ^ ~ \ bg Tr 6 ° = bg ( (Tre^)V r 2 (Tre^)V2 



<WTre A -( B+c »/ 2 < 



.4 



< log 



B + C 



Tr e J 



r£ T e B/2 e C/2 



for arbitrary Hermitian matrices A, B, C. Applying this to Equation ( I6.25P gives 

1 



9% K {n)< ^\\AH e + AH_ e \\. 
We now write out the terms in the rotated Hamiltonian in Equation f 1 6 . 2 2 j) as 



(6.26) 



(6.27) 



R e 1 H^' v R e 



Re 1 (^PAQa) + + eff x ) ■ P t (k x )\ R 6 

V A zer / 



(6.28) 



where Pt{Qa) an d Pt{k x denote the terms appearing in the Hamiltonian with periodic 
boundary conditions, obtained by mapping each site outside V to the corresponding site 
within T. 

Since k are vector operators (recall Equation f)3.12p ). 



ff x ■ (Jl e 1 n x R^j = fj x ■ R 6 (k) x = [R 1 (ff) x ] 



(6.29) 



inside T l there is no rotation, and outside we are performing an average with respect to an 
isotropic distribution, so this term makes no contribution to AH. 
As for the fixed field terms, we have 



K x . 



(6.30) 



The choices of p and 9 were intended precisely to make R x e = — e for cIl{x) > L; and for 



e • k 



the remaining (3L) d sites we have e- (^R e l K x R i 
contribution to AH which is uniformly bounded in norm by 2(3L) d 5. 



< 2, so these terms make a 
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We are left with the terms arising from the transformation of the nonrandom interaction. 
For legibility we let T = Tk for the rest of this section. For any A and any (arbitrarily 
chosen) x G A fl T, 



R- e Pr(Q A )Re = ( (g) e^' ^ e' 9 ^ Pt(Qa) (g) e - ff -* e P*- fl »)* 
\yeAnr J VzeAnr 

e -(dy-0 X )Py \ P r (Q A ) ( (^) e (*.-*.)A 



(6.31) 



vi/eAnr / \zeAnr 

(using the rotation invariance of Qa)- Expanding the exponentials, we obtain 

R-oPt{Qa)Ro = Pt{Qa) + Yl ~ °v) (PvMQa) ~ Pt{Qa)p v ) 

2/eAnr 



(6.32) 



where the estimate of the higher order terms uses 

\9 x -e y \<^ X - y ^< 7rdi ^ lA (6.33) 

and the observation that the nth order term in the expansion is potentially a sum of \A\ n 
terms, as well as ||-Ft(Qa)|| < The first order terms are odd in 9, and will cancel in 

AHg + AH^g, with the leading term being second order. What appears there is 

\2| 412 



1 (diamA) 2 |A| 2 m u \ _^rd-2 

Anry0 \ Ago 



£ (R-ePr(QA)Re - MQa)) = O ( ^ $ ~ WQ4 ) = 0(L" 

(6.34) 



where for non-finite-range interactions the last equality requires the assumption (1A.7I) . which 
is presented in Appendix [A] 

Then the right hand side of Inequality ( 16. 2 7ft is 

\\AH + AH-e\\ = 0(L d ~ 2 ) + 0(5L d ). (6.35) 

This provides only an upper bound on 5 , | / e jv/(^'), rather than a bound on its absolute 
value. However it is obvious from the definition ( I6.14p of g that g^M^f) = ~9l S m(v)j so the 
needed lower bound follows automatically. This proves (I6.14p . and thus, thorough (16.161) . 
the claimed bound ( 16.131) . □ 
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VII. STOCHASTIC LOWER BOUNDS ON THE LOCAL FREE ENERGY DIF- 
FERENCE 

To prove the Imry and Ma[l| phenomenon, Aizenman and Wehrj^] employed a somewhat 
generalized form of the central limit theorem which is suitable for the families of random 



variables presently under consideration 



We will use the following reformulation of 



that statement. Since the pro of is in the literature it will be omitted here; one incorporating 



11 



a correction due to Bovier[19|] can be found in Ref. 

Recalling that T x is the operation of translation by x G Z d and Av [-|?7a] is the conditional 
average conditioned on the values of {r] x } x& \ =: T)a, we present: 

Proposition VII. 1. Let r\ x be a collection of i.i.d. random variables (indexed by x G 7L d ) 
with Av |t7 x | 2+5 < oo for some 5 > 0, and let Gl be a family of real functions indexed by 
L G N, each with the following properties: 

1. AvG L ( V ) = 

2- |||Gx||| < e (with 1 1 | the Lipschitz seminorm defined by (I5.6p ) 

3. Gl{t]) depends on the values of rj x only for x G T L 

4. Av [G L (r])\r] A ] = G L '{T x r)\) whenever T_ X T L > = A C T L 
Then 

G l(v) / L d ^ 2 — > N(0, b) (7.1) 
in distribution as L — >• 00, for some b satisfying 

Av G\<b 2 < 2e 2 (Av \r] \) 2 . (7.2) 

Assumption H] in essence states that the different functions involved are essentially the 
same quantity at different scales, and that it is in some sense translation invariant. More 
precisely: the family forms a stationary Martingale, in a multidimensional sense of this 
condition. 

In order to use this result, we will need to establish some control over the conditions 
under which Av G\ > 0. To do this, we employ the following criterion (proven in Appendix 



III of Ref. 



4) 
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Proposition VII. 2. Let u be a Borel probability measure on R ; and let g be a continuous 
function with \ \\g\\ \ = 1 and 

g'{t)v(dt) > 0. (7.3) 



Then any of the following is a sufficient condition for J g(t) 2 v(dt) > 0: 

1. v is absolutely continuous 

2. g is differentiate with \ \\g'\\\ finite, and v has no isolated point masses. 

3. g is differentiate with \ \\g'\\\ finite, g' > 0, and v is not concentrated at a single point. 

To apply this we note that G\ is a function of one variable (rjo) which satisfies the 
conditions given above. In particular, it will be shown that the criteria 2. and sometimes 3. 



are satisfied at finite temperatures, in which case 



dGi 



dvo 

we will see that the weak FKG property implies monotonicity of G\ 



is uniformly bounded. In addition, 



VIII. CONCLUSION OF THE PROOFS OF THE MAIN RESULTS 

A. Existence of an infinite-system limit for free energy fluctuations 

Having established the properties we require of the infinite-system free energy fluctuation 
Gl, we now proceed to show that such a suitable quantity exists. This will be done in two 
overlapping cases. We will have occasion to use the following equicontinuity arguments: 

Lemma VIII. 1. Let f\j : R — > R be a family of functions labeled by i,j G N, each 
satisfying < 1 and fij{0) = 0. Then there are subsequences ik,ji such that 

/(^limlim/ y ( 2 ) (8.1) 

exists for all z G M. N . Furthermore the convergence is uniform on any compact domain 
E C R N . 

Lemma VIII. 2. Let fij : W N — )• R be as in Lemma lVIII. 1\ above, without assuming fij{0) = 
0. Then the same result holds assuming instead that there is some c such that \fij(z)\ < c < 
oo for all z G R^, i, j G N. 
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These two lemmata are proven by applying the Arzela-Ascoli theorem to obtain compact- 
ness of the space of sequences involved, then repeatedly applying the diagonal subsequence 
trick to obtain a single pair of sequences simultaneously satisfying all of the uniformity 



criteria. 



1. Finite temperature 



We now begin with the proof of Theorem IV. II in the case /3 < oo. Firstly, note that it is 
obvious that if the specified limits exist, then they will have the properties labeled CD and [2] 
in that theorem. To define the sequence used we proceed as follows. For x G Tl, let 



dGi - K lAv 



dr]x 



h-5 



then evidently 



and 



l0W^)l< e 



(8.2) 



(8.3) 



J L,K,5 



drjy 



(K x n y ) h + & (rj) - (K x K y ) n K (77) 



h-5 



(8.4) 



<2/3e 2 



For < 00, this means that for each L, <^> X LK ^ is an equicontinuous family of functions 
of the L d variables r]r L . We can apply Lemma IVIII.2I to find a decreasing sequence Si — > 
and an increasing sequence Kj — > 00 (by applying the diagonal subsequence trick, we can 
choose them to be independent of x and L) so that 



i) x L {rfi := lim lim s (rj) 

1 — >oo J — ^00 ■> 



5.5) 



exists, and by uniformity of convergence 



G L ( VL ) := lim lim ( G d £ K (r] L ) ~ Av G S £ K 



also exists with 



dG L 
dr] x 



(8.6) 



(8.7) 



for all L, rj, x E Tl- Property [3] and the relation of FKG to monotonicity then follow from 
Equation (18. 2p . and property [5] follows from Equation (18 ,4p . 
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To show property [61 we use Equations (I8.7P and (18.51) to write 



Av G[ = Av ipl = Av lim lim <f)\ K . 



(8.8) 



Recalling the definition of and using dominated convergence to exchange limits and aver- 
ages we have 



Av G[ = lim lim - Av (kq) k 

i— >oo K— >oo 2 



h+5i 



(«o) 



h-Si 
K 



(8.9) 



and since the random fields are i.i.d. we can apply the Birkhoff ergodic theorem 20J to 
replace the disorder average with a volume average, which is related to the derivative of the 
free energy by Proposition IIV.3[ giving 



Ave; 



5.10) 



2 \dh+ dh- 

Note that Equation (18.21) implies that, for any L' < L < K, 

<PI',k,5 = Av [<j>l, K , 5 \Vr L ,] , (8.H) 

and we apply the conditional form of the dominated convergence theorem to the limits used 
to define ip to obtain 

r K = ^W L \m]- (8.12) 

Since G & i K . was defined in terms of periodic boundary conditions and the limits used to 
obtain are independent of x, we also have 



and taking the limit (18.51) we have 

Together with Equations (18 . 1 2j) . this lets us obtain 

Avh& A ]=Av [iff*oT a \ Vv ] 
or in other words (in light of (18. 7p ) 



Av 



dGi 



dr]y 



'/\ 



w y GL ' {T ^ 



(8.13) 



(8.14) 



(8.15) 



(8.16) 



The uniform bounds proven above allow us to take the partial derivative with respect to r] y 
outside the conditional average [2l[], and the resulting expression can then be integrated to 
prove Property HI 
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2. Absolutely continuous random fields 



With /3 — oo, it is still obvious from Equation 15.51 that G 5 L K (0) = 0; Inequality (18. 3p 



implies 



< e, so we can apply Lemma IVIII.ll to obtain ( 18. 6p . with |||Gx||| < e. 
AyGl = is obvious. We can apply the diagonal subsequence trick to obtain sequences 
independent of L, which implies that the consistency condition @] of Theorem IVII.ll is satis- 
fied. 

The difference arises in proving that Av G[ = eM. Without equicontinuity of the deriva- 
tives, there is no reason to expect that functions constructed like <j) of the previous section 
will converge uniformly, or that the pointwise limit will be different iable. This problem can 
be resolved applying the following convergence criteria: 

Lemma VIII. 3. Let g be a measurable function and g n a sequence of measurable functions 
such that Halloo < 1; IMIoo < 1? an d 

r-6 r b 



lim / g n (t)dt = / g(t)dt 

n ^°°Ja J a 



(8.17) 



for all a,b £ K. Then for any finite measure v which is absolutely continuous with respect 
to the Lebesgue measure, 

lim / g n dv = g du. (8.18) 



This can be proven by elementary measure theory techniques 11]. 



Thanks to Rademacher's theorem 



22], Lipschitz continuity guarantees that f n have dis- 



tributional derivatives, i.e. functions g n satisfying 

g n {t)dt = f n {b) - f n (a) 
for any a, b £ R; and furthermore that ||g n ||oo = |||/n||| < 1- We therefore have 



(8.19) 



Corollary VIII. 4. Let f n be a sequence of functions M — > M such that f n —>f pointwise, 
an d |||/n|| I — 1- Then f has a distributional derivative f satisfying 



lim / g n (t)du(t) = [ f'(t)du(t). 
Applying Corollary IVIII.4I twice to G^ k .(vl) — AvG^. and G x , we obtain 



(8.20) 



AvG; = lim lim Av [G? K . (rj) 



(8.21) 
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and applying Proposition IIV.3I and the Birkhoff ergodic theorem to the right hand side as 
in the finite temperature case we obtain 

AvG-^-iq (8,2) 



2 V dh+ dh- 
and the proof of Theorem IV. II is complete. 



B. Concluding the proofs of Theorems IIV.1I and IIV.2I 

In summary, we have established that under the hypotheses of Theorems II V. ll and II V. 21 1 he 
free energy fluctuations Gi,(rj) are uniformly bounded above by CL d ~ 2 , and correspondingly 
CL d / 2 , and that if there is a 1 st order phase transition (i.e. a discontinuity in then Gi{rj) 
has asymptotically normal fluctuations on a scale of L d / 2 . That is is a contradiction even 
in the marginal dimensions (where the powers match) as a normally distributed random 
variable will take with positive probability values which are arbitrarily large on any given 
scale. An alternative way to present this case is to note that assuming the existence of a 1 st 
order phase transition, by Proposition IVII.ll we get 

lim Av exp (tG L /L d/2 ) = exp{t 2 b 2 /2) . (8.23) 

At the same time, if ||Cr.r,||oo < AL d / 2 then for all positive t: 

A ve tG ^ Ld/2 <e tA . (8.24) 

Clearly if b ^ 0, the two relations are incompatible at sufficiently large t. Finally Propo- 
sition IVII.2I allows to conclude that in any of the cases under consideration, 6 = implies 
M = 0, and the proof is complete. 



Appendix A: Decay conditions for long range interactions 



In this Appendix we turn to the general class of interactions under which the results 
derived in this work apply. 

The estimate which is relevant for Theorem IIV.1I is valid if: 



diamA^||Q A || <dL [2 ~ d ^ 2 (A.l) 
Abo ' ' 

diam A<L 
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for some constant c'. In case of two body interaction, that is when Qa is non zero only when 
A is a two point set, this condition is met if 

\\Q{x,v}\\ < Const./lx-yl 3 ^ 2 . (A.2) 

As noted in the proof, the essential estimate which is required for Lemma IVI.2[ and 
therefore Theorem IIV.lt is that 

Av\\V^ L \\<cL d - 1 + bL d/2 , (A.3) 

with some finite constants b and c. The contribution in VJl of terms involving the random 
fields 77 consists of a collection of OIL' 1 " 1 ) independent summands with finite average norm. 



Thus the assumption (IA.3[) requires only a bound o 



of the interaction. The following simple calculation 
condition for that. 

By the triangle inequality 



: the same form on the nonrandom part 
4j shows that flA.lj) provides a sufficient 



\K'°\\< (A- 4 ) 

In this sum the terms with diameter L or less contribute at most 

2dL d - 1 ^^<2dc'L d /\ (A.5) 



A30 
diam A<L 



and the remaining portion is bounded by 



A30 
diam A<L 



<L d - 1 ^ diamAi^lig^H <c'L d ' 2 

A30 ' ' 

diam A<L 



(A.6) 



Putting the two parts back together we have Inequalit y (A.3P with C = (2d + l)c' '. 



The recent result of Cassandro, Orlandi and Picco [23[ allows to conclude that at d — 1 
the condition ( 1A.2I) does indeed provide the threshold decay rate for the validity of Theo- 
rem IIV.1I They prove that the phase transition is stable under weak disorder in a family of 
one dimensional Ising spin systems with long range interactions with decay rates arbitrarily 
close to 3/2, more explicitly with ||Q{ x ,j(}|| ~ Const. j\x — y| 3 ^ 2-7 at any 7 £ (0,0.08). 
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(Curiously, while the proof of the complimentary result does not extent beyond this narrow 
range of values of 7, the lower end of the interval coincides with flA.2j) .) The question 
whether the above criterion is optimal also for d — 2 is of interest, as there are some notable 



systems with inverse cube interactions in two dimensions 



24]. 



For the stronger conclusion which is derived here for interactions with continuous sym- 
metries the relevant assumption is not (1A.3|) but 



^(diamA) 2 |A|IIQA|| < 00, (A.7) 

which is used in (16.341) , in the proof of the free energy estimate of Lemma IVI.3I For pair 
interactions this reduces to the statement 

\\Q{o,x}\\ M 2 < 00 > ( A - 8 ) 

found (in slightly different notation) in Ref. If]. 
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